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ABSTRACT 

O . 

O . The rapid variability of the VHE emission reported for some TeV blazars 

^ ! implies Doppler factors well in excess of those inferred from superluminal motions 

^! and unification schemes. We propose that those extreme flares may result from 

^ ' radiative deceleration of blobs on scales where local dissipation occurs. The 

minimum jet power estimated from the resolved synchrotron emission on VLBI 
scales appears to be consistent with this model. It is shown that if the energy 
^ i distribution of nonthermal electrons accelerated locally in the blob is reasonably 

flat, then a background radiation field having a luminosity in the range 10'^^- 
2 ! 10"^^ erg s^^ can give rise to a substantial deceleration of the blob, but still be 

to ! transparent enough to allow the TeV 7-rays thereby produced to escape the 

system. 

^ i Subject headings: galaxies: active - quasars: general - radiation mechanism: 

T^ij- i nonthermal - X-rays:galaxies 

in 

O I 1. Introduction 

o 



X 



Very high energy emission (VHE , -E ^ 100 Ge V) has been detected from over a dozen 



blazars (for an updated list see, e.g.. IWagner 20071 ). all of which are exclusively associated 



^ I with the class of high peak BL Lac objects. The observed bolometric luminosity during 
quiescent states is typically of the order of a few times lO"^"^ ergs s~^, with about 10 percents 
emitted as VHE 7-rays. The luminosity in the VHE band may be larger by a factor of 10 to 
100 during flaring states. The intrinsic spectra (corrected for absorption on the extragalactic 
background light) appear to be rather steep, but seem to harden with increasing flux. A peak 
photon energy in excess of 10 TeV has been measured in the most extreme cases. Dev iations 



from a power law are clearly seen in the data in some cases (e.g.. lAlbert et al. 20071 ). 



Rapid variability is a characteristic property of VHE blazars. Large amplitude vari- 
ations of the VHE 7-ray flux on time scales of several hours and less have been reported 
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for Mrk 421, Mrk 501 and PKS 2155-304. This rapid variability implies that the emis- 
sion at the observed energy originated from a region of size (as measured in the Lab 
frame) Ar ^ 10^'^rDtvar,h/(l + z) cm, where F and D are the bulk Lorentz factor and 
the corresponding Doppler factor of the emitting matter, respectively, tvar,h is the observed 
variability time in hours, and z is the redshift of the source. If the emission originated 
from small jet radii, rem ~ Ar, then the requirement that the 7 rays will not be ab- 
sorbed by pair production on local synchrotron photons implies hig h Doppler factors, D ~ 
30 — 100 (IKrawezynski 20021 : iLevinson 20061 : iBegelman et al. 20071 ). Such high values are 
consistent with those obtained from fits of the SED to a homogeneous SSC mo but 
are in clear disagreement with the much low'e r values inferred from unification schemes 
( Urrv fc Padovani 199ll: iHardcastle et al. 20031) and supe rluminal motions on parsec scales 
f lMarscher 19991 : iJorstad et al. 200ll : iGiroletti et al. 2004J'). Various explana tions, including 
a struct ure consisting of interactin g spine and sheath ( Ghisellini et al. 20051). opening angle 



effects ([Gopal-Krishna et al. 20041 ) and jet deceleration (iGeorganopoulos fc Kazanaz 2003 



Finer fc Edwards 20051 ) have been proposed in order to resolve this discrepancy. A variant 
of the decelerating jet model is considered further below. 

It has been argued that such high values of the Doppler factor may not be required if 
the 7-ray production zone is located far from the black hole, at radii rem >> Ar. In that 
case the compactness of the TeV emission zone may be constrained by the variability of the 
IR flux observed simultaneously with the TeV flare, allowing low values of D in cases where 
the variability time of the IR emission is much longer than the duration of the TeV flare. 
However, the fraction of jet energy that can be tapped for production of 7 rays in a region 
of size Ar located at a radius rem is 77 ~ {Ar/9rem)^- As a consequence, either the opening 
angle of the jet must be very small, 6 ~ Ar/rem << 1, or the jet power must be much larger 
than the luminosity of the TeV emission measured during the flare. Another possibility 
is that the TeV emission is produced by a converging sh ock in a reconflnement nozzle, as 
proposed for the HSTl knot in M87 (ICheung et al. 20071 ). It shoul d be noted though that 
in M87 the X-ray a nd TeV luminosities, Ltby ~ -f/x ~ 10"^^ erg s~^ f lAharonian et al. 2006 



Cheung et al. 20071 ). are much smaller than the TeV luminosity, LTev ~ lO^'^ erg s ^, 



observed typically in the TeV blazars. Estimates of the je t power in M87 yield Lj ^ 10 erg 



s [e.g. 
fraction, 



Bicknell fc Begelman 19961 : IStawarz et al. 20061 ) , implying a very small conversion 



L^ev/Lj ^ 10 . Even with such a small conversion efficiency an opening angle 
6 < 10~^ rad is required if the TeV emission were to originate from the HSTl knot, unless 
reconflnement c an give rise to suffic ient convergence of the jet at the location of HSTl, 
as proposed by (ICheung et al. 20071 ). We note that even modest radiative coohng of the 
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shocked jet layer in a proton dominated jet will lead to such a convergence, at least in the non- 



relativistic case (jPeter &: Eichler 19951 ). The effect of cooling on the coUimation of relativistic 



jets needs to be explored, b ut we expect a sirn ilar behavior. Stationary radio features 



observed on VLBI scales (e.g.. |Jorstad et al. 200ll ) seem to indicate that recollimation shocks 



may be an i mportant dissipation channel in blazsrs, and this may apply also to other sources. 



e.g., GRBs (IBromberg fc Levinson 20071 ). Whether the extreme TeV flares observed in VHE 
blazars can be accounted for by recollimation shocks at radii rem >> remains to be 
explored. As stated above, this would not resolve the 'Doppler factor crises' if the IR 
emission will turn out to vary on timescales comparable to the TeV emission. 

( Georganopoulos fc Kazanaz 20031 . GK03) proposed a scenario in which the deceleration 



of the fast jet base is mediated by Compton scattering of synchrotron photons produced 
further downstream, in the slow part of the flow. The major fraction of the bulk energy 
is assumed to be converted to TeV radiation within the transition layer, so there is no 
missing energy problem in this model. The main motivation in that paper was to reproduce 
the observed SED in a source that propagates at mildly relativistic speeds on VLBI scales. 
Using a given bulk velocity profile for the decelerating plasma GK03 computed the spectrum 
emitted from the jet and argued that deceleration from a modest Lorentz factor (F ~ 15) 
down to r ~ a few in a jet observed at sufficiently small viewing angles (6'„ ^ 3°) can 
indeed account for the average SED observed in TeV blazers. The dynamics of the system 
has not been treated in a self-consistently manner in GK03. In particular, it has not been 
demonstrated that (i) the backward emission from the slow jet section can provide sufficient 
radiative drag to decelerate the fast jet, and (ii) that the TeV photons can escape to infinity. 
Moreover, as explained above the Doppler factors inferred from opacity constraints during 
the extreme flaring states observed in Mrk 412, Mrk 501 and Pks 2155-304 are considerably 
larger than those invoked in GK03 to explain the average broad-band spectrum of these 
sources. 

In this paper we consider the possibility that those extreme flares are produced by radia- 
tive deceleration of fluid shells on scales where local dissipation occurs (r^ ~ 10^ — lO^r^). The 
dissipation may be accomplished through formation of internal shocks in a hydrodynamic jet 
or dissipation of magnetic energy in a Pointing flux dominated jet. A similar model has bee n 



proposed earlier for flares in EGRET blazars (IRomanova fc Lovelace 19921 : iLevinson 19981 ). 



For the TeV blazars a background luminosity Ls ~ 10'*^ — 10'*^ erg s~^ would lead to a 
substantial deceleration of the front and still be transparent enough to allow the TeV 7-rays 
produced by Compton scattering of the background photons to escape the system, provided 
the energy distribution of radiating electrons is sufficiently flat. The ambient radiation field 
is most likely associated with the nuclear continuum source. The bulk Lorentz factor of 
the jet during states of low activity may be appreciably smaller than that of fronts expelled 
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during violent ejection episodes. 



2. Minimum Jet Power: A Consistency Check 

The decelerating jet scenario implies that the major fraction of the bulk energy is radi- 
ated away in the form of VHE photons on very small scales (~ 10^ — lO^r^). Consequently, 
the remaining jet power on VLBI scales must be much smaller than the luminosity of the 
VHE 7-ray emission. As a consistency check, we estimate the minimum power of the pc 
scale radio jet. From the resolved radio synchrotron emission emitted at a radius r = Irpc pc 
from a jet of opening angle 6 = 0.1^_i rad we define a fiducial equipartition magnetic field 
strength B* ~ b{TBg/ O^iTpcY^'^ mG, where the brightness temperature Tb = lO^Teg K 
is e valuated at the observed freq uency z/ = z/g GHz, and the radio spectral index is ~ 0.5 



(cf., I Rybicki fc Lightman 19791 ). Let the jet Lorentz factor be T and let its velocity makes 
an angle 9n with the line of sight. The jet power associated with the emitting electrons and 
the electromagnetic field is then given by 

c [TB*9r\^ 40 / 6'-iV^^°^^r2T^4/7, 5/7 



^^•^il'^J ~10-^:^J r^T^^r ergs- (1) 
where D = [T{l—j3 cos On)]~^ is the Doppler factor and rough equality occurs at equipartition. 



VLBA images of Mrk 421 at 15 and 22 GHz reveal a core jet morphology (iMarscher 19991 ) 
with some jet components at a linear distance of Vpc ~ 1 from the core. The jet components 
contain only a few percents of the total flux density (of the order of 1 Jy). The resolu- 
tion at 15 GHz is about 0.5 mas. From the above we estimate Tb9(z/ = 15GHz) ^ 1 for 
the jet, and a minimum power of -Lj,min ^ 10^^ erg s~^. For Mrk 501 we obtain similar 
numbers. Thus, the minimum jet power inferred from the resolved radio emission is well 
below the TeV luminosity measured in these sources, L^ev ~ 10 ^4 qj.^ g-i_ "Yhis should be 



contrasted with FRII sources where Tg ^ 10^^ K are measured (IReadhead 19941 ). implying 
a minimum jet power close to the Eddington limit, and radio jets in microquasars, where a 
minimum jet power in excess of the Eddington limi t has been measured on scales of ~ lO^r^ 



( ILevinson fc Blandford 19961 : IPistefano et al. 20021 ). 



Radiative Deceleration of Relativistic Shells 



Consider a fluid shell, described by a stress-energy tensor 



T"^ = hnu^'u^ - pr]""^ + ^ + ^r^^^F^) , (2) 
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where F^^j is the electromagnetic tensor, = (F, T(3) is the 4-velocity of the bulk fluid, 
and ra, p, and h, are the proper particle density, pressure and specific enthalpy, respectively, 
interacting with some ambient radiation field. Suppose that the shell has been ejected and 
accelerated to a bulk Lorentz factor Fq at some radius ~ lO^r^, at which dissipation 
suddenly commences, e.g., due to collision with another shell or with a confining medium. 
The dynamics of the front is governed by the energy momentum equations. 



d 



T 



afj, 



Si;, 



(3) 



where the source term Sj^ accounts for the radiative force acting on the front. In terms of 
the distribution functions of the target photons, fs, and electrons (we don't distinguish here 
between electrons and positrons), /e, the source term associated with the radiative force is 



given, in the limit of Thomson scattering, by (jSikora fc Wilson 198ll : iPhinney 19821 ) 



= -car 



Pi 



P'i 



' fsfePeaPs 



P^ + 



{PeuPs)Pe 



(4) 



with and being the 4 momenta of electrons and soft photons, respectively, as measured 
in the Lab frame. We suppose that the electron distribution function is isotropic in the 
fiuid rest frame and can be approximated as a power law: dnJ^jd'y = ATcnieCp'^ fe{p) oc 7"''; 
7mm < 7 < 'Jmax, whcrc uieC^'^ is the Corresponding electron energy, as measured in the 
comoving frame. We further suppose that the photon distribution is isotropic in the star 
frame. Under the above assumptions the zeroth component of eq. (jlj) yields. 



Here 



fed'pe = I ^d7, 



is the proper number density of nonthermal electrons, 

< 7' >= — / l^fed^Pe, 



(5) 
(6) 

(7) 



and Us{r) is the total energy density of the target radiation field at radius r. Let u'^ = 
J '-)mf.(? fed^Pe be the proper energy density of nonthermal electrons, < 7 > mgC^ = Vtlju^ 
their average energy, and define < 7^ > / < 7 >= X7max- In terms of these quantities the 
radiative force term (eq. [5]) can be re-expressed as 



5° 



4(T7 



Xr'7 



3 ' 
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For the power law energy distribution invoked above we have x = (2 — q')/(3 — g) if g < 2, 
X = [Mlmax/lmin)]''^ if ? = 2, and X - {iminhmaxY''^ if g > 2. With 7„j„ = mp/me and 
Imax = 10^ we have x > 0.1 for g < 2. 

The energy flux of the decelerating front can be obtained from eq. (j2]): 

T^*^ = {nh + fi'V47r)r2/3 = ^T^/?, (9) 

where n and h are defined above and B' is the rest frame magnetic field. Using eqs. ([3]), ([5]) 
and (IHD we arrive at, 



For illustration suppose that Us(r) oc and that the proper density and average energy of 
the nonthermal electrons are independent of radius. In the limit /3 = 1 the solution to eq. 
f fTU]) reads 

- - ' (U) 



i + 

where Fq = r(r = r^), and the stopping length / is given by 

3meC^ , , 

I = — Tj^ — r^- (12) 

In the last equation = u'^/u'j < 1 denotes the fraction of total jet energy carried by the 
nonthermal electron population. 

In order to calculate the pair production opacity the spectrum of the target radiation 
field must be specified. For simplicity let us assume a power law spectrum, Isi^s) oc e^"; 
^s,min < Cs < es,max, whcrc Is = hcp^fs is the intensity and = ps/meC denotes the photon 
energy in rrieC^ units. The optical depth for absorption of 7 rays of dimensionless energy 
by pair creation is then given to a good approximation by 

/meC^)e^g{e^), (13) 

with g{e^) = (e^es,™„)""^ if a > 1 and g{e^) = {e^es,maxT'^ if a < 1, and g{ej) < 1 in both 
cases. Using eqs. (fT2l) and (fT3|) . the stopping length can be expressed in terms of the 7-ray 
optical depth as, 

1 = 1 (^) ( ^) (14) 

rd XG^77 \(^tJ Vro7max/ 

Note that since the observed 7-rays are produced by Compton scattering of the front electrons 
we must have Fo7max > £7 for any e^. The highest 7-ray energy observed is likely to be limited 
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by opacity. Denoting by e^^th the photon energy at which the pair production opacity is unity, 
viz., r^^(e^^th) — 1, and taking a^^/ax = 0.2 we find that //r^^ < 1 if 

Adopting for illustration x = = 0.1, g{e^^ih) = 0.1, we conclude that extension of the 
nonthermal electron spectrum to a maximum energy 7max of the order of a few times e^^th/ro 
is sufficient to cause appreciable deceleration of the front. Gamma rays having energies 
above e^^th will be degraded to somewhat lower energies by virtue of the relatively large 
pair production opacity. Thus, we naively anticipate some accumulation of fiux around 
the maximum observed 7-ray energy, at ~ a few TeV, although we stress that detailed 
calculations are required to determine the exact shape of the spectrum. The corresponding 
luminosity of the ambient radiation can be estimated from eq. f fT3|) to be 

Ls ^ "^'"'Y ^ ^ 10^%i6(6,,th/10^)^?-^ erg (16) 

Assuming the acceleration rate of nonthermal electrons to be on the order of their gyro- 
frequency, we estimate 7max ~ 10'''^ for magnetic field energy density u'^ ^ O.lu^-. For Fq = 30 
this corresponds to Fo7max/e7 — 50 for the highest 7-ray energy observed (~ 10 TeV). If the 
ambient photons originate from a central continuum source (e.g., accretion disk), and about 
10 percents are intercepted by the jet, then of the order of 10^^ — 10^^ erg s~^, roughly the 
luminosity observed in LLAGNs, is required to accommodate the spectrum and variability 
of the TeV blazars. 

We consider now the possibility that the target photons originate from the slow part 
of the flow (GK03). Assuming for illustration that the proper energy density u'j remains 
constant we have Ljoo = (Foo/Fo)^-Z^j. The fraction of synchrotron luminosity emitted back- 
wards (into a solid angle in the rest frame of the slow flow), as measured in the Lab 
frame, is ~ F^"^. Consequently, the luminosity of target photons is at most V^Lj^ or 
Ls ^ Lj/(FooFo)^ Adopting Lj = 10^^ erg s ^, roughly the observed TeV luminosity, and 
Foo = 4, we deduce that Fq cannot be larger than 20 or so in order that Lg will satisfy 
the requirements imposed in eq. f|T6|) . Note also that this photon mediated transition is 
unstable, so that some additional mechanism is required to maintain the flow at low T^o 
downstream. 



4. Discussion 

The rapid variability of the VHE 7-ray emission reported for some VHE blazars implies 
Doppler factors D 30, much larger than those inferred from superluminal motions and 
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unification schemes. A plausible mechanism that can generate such flares is radiative decel- 
eration of high Lorentz factor shells. The minimum jet power estimated from the resolved 
synchrotron emission on VLBI scales is consistent with this model. Deceleration to F ~ a few 
can be accomplished if a signiflcant fraction of the shell's bulk energy is carried in the form 
of nonthermal electrons with a sufficiently flat energy distribution [dnejd'y oc 7"'^; g < 2) 
that extends up to a maximum energy at which the pair production opacity is of the order 
of a few or larger. The spectrum emitted during the flare is expected to be hard, with vFy 
peaking roughly at an energy at which the pair production opacity is unity. Deviation from 
a power at the highest energies are expected, owing to attenuation by pair production an d 



temporal effects. Spectral curvature is indeed observed in some cases (lAlbert et al. 20071 ). 



Correlation between the VHE 7-ray emission and emission in other bands, particularly X-ray, 
is naively anticipated, although detailed calculations are required to assess exact relations 
(e.g., time lags and amplitude ratios). Such correlations have been reported for Mrk 421 



(IFossati et al. 20041 ). The properties of VHE flares produced by this mechanism should dif- 



fer from those predicted for the pair cascade jets in p owerful blazars, like 3C279, where 



propagation from low-to-high 7-ray energies is expected (IBlandford fc Levinson 19951 ). 



I thank C. Dermer for useful comments. This research was supported by an ISF grant 
for a Israeli Center for High Energy Astrophysics. 
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